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Abstract 

We propose a curvaton model in which the initial condition of the curvaton 
oscillation is determined by its attractor behavior during inflation. Assuming a 
chaotic inflation model, we find that the initial condition determined by the attractor 
behavior is appropriate to generate a sizable non-Gaussianity contribution to the 
curvature perturbation, which will be tested in the foreseeable future. Implications 
on the thermal history of the universe and on particle physics models are also 
discussed. 



1 Introduction 



The origin of the large scale structures of the universe and the fluctuation of the cosmic 
microwave background (CMB) radiation can be successfully explained by the primor- 
dial density perturbation generated by cosmic inflation [1] [2] . In particular, the simplest 
scenarios in which the density perturbation is dominated by the fluctuation of a single 
inflaton field has been very successful so far. 

In general, however, any light fields other than the inflaton also fluctuate during 
inflation and may contribute to the primordial density perturbation. In fact, in the 
curvaton mechanism j3] , the primordial density perturbation is mainly generated not by 
the inflaton fluctuation but by the fluctuation of a light field, the curvaton. The curvaton 
model is attractive since the inflaton field does not need to provide the whole density 
perturbation any more, which relaxes the constraints on the inflation models. 

One of the interesting observable feature of the curvaton scenario is that it can lead 
to large non-Gaussianity in the density perturbation when the curvaton energy density is 
subdominant at its decay time. This is in contrast to the primordial density perturbation 
generated by the inflaton fluctuation, where the non-Gausssianity is predicted to be highly 
suppressed [3] . Therefore, if the large non-Gaussianity suggested by WMAP data j5j E] 
is confirmed by the forthcoming Planck experiment [7], the curvaton scenario becomes a 
more plausible explanation for the density perturbation. 

There is, however, a drawback in the curvaton scenario if it generates sizable non- 
Gaussianity. That is the initial condition problem of the oscillation of the curvaton field. 
As we will briefly review in the next section, the initial field value of the curvaton oscil- 
lation is required to be about a hundred times of the Hubble parameter during inflation 
to generate sizable non-Gaussianity. This initial field value of the curvaton is, however, 
difficult to be justified since there seems no special meaning on such an intermediate 
scale E] 

1 If the required initial field value were at the origin of the curvaton field or at around the Planck scale, 
for example, it could be rationalized. That is, the origin of the curvaton can be a symmetry enhancement 
point, which can be chosen during inflation if the symmetry is also respected by the inflaton dynamics. 
The field value at the Planck scale can be also chosen if the field value during inflation is determined 
stochastically [9] , since the scalar potential of the curvaton is expected to become very steep for the field 
value larger than the Planck scale. 
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In this paper, we propose a curvaton model in which the initial condition of the 
curvaton oscillation is set by the so-called attractor behavior (See Ref. [8], for example) 
during inflation. As we will show, the initial condition fixed by the attractor behavior can 
be appropriate for the curvaton scenario to account for the non-Gaussianity suggested by 
the WMAP data. We also discuss how high the decay temperature of the curvaton field 
can be. The implications on the particle physics models are also discussed briefly. 

This paper is organized as follows. In section[2j we briefly review the curvaton scenario. 
In section[3j we discuss the curvaton model which has the attractor behavior. There, we 
also discuss some implications of the model on the particle physics models. 

2 Brief Review of The Curvaton Scenario 

Before going to discuss the attractor behavior of the curvaton field, let us briefly review 
the basics of the curvaton scenario. For simplicity, we assume a simple quadratic potential 
for the curvaton field, a, 

V(a) = \ml^ , (1) 

where m denotes the mass of the curvaton. In the curvaton scenario, the curvaton mass is 
assumed to be smaller than the Hubble parameter during inflation. We also assume that 
the Hubble induced mass is suppressed. With such a flat potential, the curvaton field is 
over-dumped and its field value is frozen to its initial value, <7j, during inflation. It should 
be noted that the field perturbation around <jj, o = a; t + 5cr, becomes the origin of the 
primordial density perturbation of the universe eventually. 

After inflation, the curvaton starts to oscillate around its vacuum value (cr = 0) once 
the Hubble parameter falls below the mass of the curvaton, i.e. H < m a . Eventually, 
the coherent oscillation of the curvaton decays into radiation at t = tdec- Hereafter, 
we assume that the reheating process after inflation completes before the decay of the 
coherent oscillation of the curvaton. 

In the curvaton scenario, the curvature perturbation ( on spatial slices of uniform 
density is evolving until the decay of the curvaton (see e.g. Ref. [H] for a review). After the 
decay of the curvaton oscillation, the curvature perturbation stops evolving and becomes 
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constant at the super- horizon scale, which is given by [TU] . 
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where p r denotes the energy density of the radiation, p a and 5p a the energy density of the 
curvaton and its fluctuation on the spatially flat slice. In the above expression, rd ec is the 
ratio of the energy densities at the decay time of the curvaton, r dec = Paitdec)/ Pritdec), and 
we have used the fact that Sp a /p a is constant in time from the horizon exit to the decay 
time. It should be noted that the inflaton contribution to the curvature perturbation is 
assumed to be negligible. 

By using Eq. (pi), we obtain the power spectrum of the curvature perturbation, 



1 dec 

16vr 2 



Hk 



(3) 



(4 + 3r dcc ) 2 \27rcTj 

for rdec *C 1. Here, Hk is the Hubble parameter at the horizon exit of the wave number k 
and we have used Sa^ = Hk/2-K. The nonlinearity parameter f^i is, on the other hand, 
defined by the bispectrum, 

6 



£ c (ki,k 2 ,k 3 ) = -fNL(ki,k 2 ,h)\P c (ki)T> c> (k2) + cyclic permutations] 



(4) 



which can be extracted from Eq. ^ by remembering the relation to the Gaussian pertur- 
bation (g, 



(5) 



(6) 



C - Cg + -^fNLCg ■ 

By comparing Eqs. ^ and ([5]), we find that 

>. r dec H k 5 

Cg - ~ A , /nl ^ , , 

47T o-j 3r dcc 

for r de c <C 1- The above expression shows that the nonlinearity parameter f^L can be 
sizable if the curvaton energy density at its decay time is subdominant. 

Now, let us estimate the ratio of the energy densities, r dec . Here, we assume that the 
curvaton starts to oscillate before the reheating. At the beginning of the oscillation [H] 

set by, 



H 2 ~ 



1 dV 
5a da 



Hi 



(7) 
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the ratio of the energy densities of the curvaton and the inflaton, 0, is given by, 

P<t> SM^L 6 Ml ' 1 ] 

which is a constant of time until the inflaton decays and the universe is reheated. Here, 
Mpi denotes the reduced Planck scale. Once the reheating process completes, the energy 
density of the inflaton is converted to the radiation energy density which is diluted faster 
than the curvaton density. As a result, the ratio of the energy densities at the curvaton 
decay time is given by 

5 erf T R 



6 Ml T dec 



(for T R > T dec ) , (9) 



V 

where Tr and Td ec are the reheating temperature and the decay temperature of the cur- 
vaton, respectively. 

As a result, we find that the initial value of the curvaton oscillation should be much 
smaller than the Planck scale, i.e. a, <C M p i, for a sizable nonlinearity parameter to be 
generated. More explicitly, we may express the required initial value <7j in terms of the 
power spectrum and the nonlinearity parameter by using Eqs. Q and (§), 

5 H k 




(10) 

where k* = 0.002Mpc _1 is the pivot scale. Here, we have again assumed r^ec *C 1. 
As mentioned in the introduction, this initial value of the curvaton is difficult to be 
rationalized, since it seems no apparent physical meaning on that scale. As we will 
discuss in the next section, however, the above initial condition is dynamically realized in 
the model with the attractor behavior. 

Before closing this section, let us comment on the decay temperature, T dec . First, let 
us remember that the baryon asymmetry of the universe and dark matter are required to 
be generated after the curvaton decay to avoid unacceptably large iso-curvature pertur- 
bations [TU]. Thus, if the baryon asymmetry is explained by thermal leptogenesis [TS], for 
example, the decay temperature of the curvaton is required to be high, Td ec > 10 9 GeV, 
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since thermal leptogenesis requires a high temperature environment O(10 9 ) GeV [16J. On 
the other hand, from Eqs. (|6]), (|9| and (10), the decay temperature is required to be 

25 1 / H k V 
dec " 28^V c {k)f NL \M^J *' 



5xlO™GeVJ \f NL J \vl /2 {K) / Vl0 13 GeV 



c 

Thus, a rather high decay temperature, Td cc > 10 9 GeV, for example, can be achieved 
if the Hubble scale during inflation and the reheating temperature are rather high for 
f NL = O(10). 

Unfortunately, however, these requirements are not easily satisfied in explicit inflation 
models. To make our discussion concrete, let us consider the chaotic inflation model [TT] 
with a quadratic potential ^m|0 2 as an example. In this model, the Hubble parameter is 
given by = \j2Nk/3m ( f ) where is the number of e-foldings, and hence, the above 
Tdec is reduced to 

rp 25 N k ' 2 

-L dec 




432ir*V c (k)f N L 

It should be noted that the inflaton mass cannot be much larger than = 10 12 GeV, 
since otherwise the inflaton contribution to the curvature perturbation, 

N 2 rn? 

W*) mflat0n - , (13) 

cannot be ignored and the non-Gaussianity is suppressed]^] Therefore, the decay temper- 
ature is expected to be rather low for a sizable nonlinearity parameter, /jvl = O(10), 
unless the reheating temperature after inflation is much higher than the inflaton mass 
(see Refs. [19] for discussions on the reheating temperature and the inflaton mass). 



2 If the inflaton contribution to the curvature perturbation is sizable, baryogenesis before the curvaton 
decay can be consistent with the current constraint on the isocurvature fluctuation |18j . 
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3 Curvaton Model with Attractor Behaviors 



Now, let us discuss the curvaton model in which the curvaton field value after inflation 
is set by the attractor behavior. To be concrete, we consider the chaotic inflation model 
with a curvaton field whose scalar potential is given by, 

V{M = \m%<p + + \mla 2 . (14) 

where A denotes a coupling constant. We assume that the Hubble induced mass is neg- 
ligible. Notice that the spectral index n s of the curvature perturbation is red-tilted 
(n s ~ 0.98) in the chaotic inflation scenario, which is consistent with the CMB obser- 
vations [H]. 

3.1 Dynamics of the fields during inflation 

When the slow-roll conditions are satisfied, the dynamics of the inflaton field and the 
curvaton field are described by 

3H<j) = -d^V, 3Ha = -d a V . (15) 

In the followings, we assume that the initial value of the inflaton at the beginning of 
inflation, <fio, is much larger than the Planck scale since we are considering the chaotic 
inflation scenario. The initial value of the curvaton, <7o, is, on the other hand, at around 
the Planck scale, which can be justified if the curvaton potential steeply increases for 
° > M pl . 

In the large curvaton field value region, a > \^2m (7 /\/X = cr 2) the potential of the 
curvaton field is dominated by the quartic term. In this region, the solutions of the 



equations of motions of the inflaton and the curvaton (Eq. (15)) are interrelated by, 

1 . 4) = J_ ( l_ _ 1 



2 log7 = ^h" , (16) 
mi 2A \ cr (Xn 



where the subscript "0" denotes the initial values at the beginning of inflation. After 
enough time, a becomes far smaller than <7o, and hence, the curvaton field enters into the 
attractor solution given by 

-1/2 

X 2A 



= ^ log? ) • (17) 
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Figure 1: The attractor behavior of the curvaton field during inflation. The magenta, blue, 
green and red lines correspond to the initial condition <jq = 1, 10 _1 , 10~ 2 , 10~ 3 xM p i, respectively, 
while the initial condition of <po is 10 2 , 10 3 , 10 4 x M p \ from left to right. The figure shows that the 
curvaton becomes insensitive to its initial condition. In the colored region, slow-roll condition 
for the curvaton held is violated and the curvaton field promptly rolls down until it gets out of 
the colored region, as is shown with the dotted line. Afterwards, the curvaton field follows the 
attractor solutions shown in the figure. 

In Fig.[IJ we show an attractive behavior of the massless curvaton during inflation for 
a small initial value o"o < M p \. The figure shows that the curvaton evolves along the 
attractor solution and becomes insensitive to its initial value o~o- 

When the curvaton field becomes smaller than a < a 2 during inflation, the evolution 
of the curvaton deviates from the attractor solution, while it keeps following the attractor 
behavior by the end of inflation if the curvaton field satisfies o > 02 until then. In both 
cases, the initial condition of the curvaton oscillation after inflation, <Tj, is insensitive 
to its initial condition at the beginning of inflation, o~o- As we will show shortly, the 
dynamically chosen o~i can be appropriate to generate a large /nl, and hence, the model 
with the attractor behavior provides us with a solution to the initial condition problem. 

Before proceeding to the details of the curvaton behaviors, we comment on the initial 
condition of the inflaton field. It should be noted that the initial condition of the inflaton 
cannot be arbitrarily large since the stochastic behavior is significant for (pH^ 1 < H/2tc. 
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For the quadratic potential, this condition corresponds to 

> (967ry/ 4 M p 3 / 2 K /2 ~ 10 4 x M pl / ^ G, " n * 



J 

If the universe begins when the inflaton escapes from the stochastic region, the initial 



condition is given by the field value which saturates the inequality of Eq. (18). In the fol- 
lowings, we assume that the initial condition of the inflaton field is given by the saturated 
value, 

0o ~ 10 4 x M pl f^H^L) 2 ( (19) 

although the attractive behavior is less sensitive to the initial condition of the inflaton, 
0o (see Fig. [I). 

3.2 cr 2 < cr att (0 C nd) 

When the curvaton mass is small enough, it follows the attractor solution by the end of 
inflation, 

~ 2M pl = end . (20) 
Thus, in this case, the initial condition of the curvaton after inflation is given byj^] 

<Ti ~ (T cnd = -j== log" . 21 

V2A V 0cnd/ 

After the end of inflation, the inflaton field starts oscillation around the minimum of the 
potential. By remembering that 



2 11-22 ^ 2 



1 dV 
5cr da 



g^end — ^Q m 1 X end < ^end 5 (22) 



A Q 1 

end 

at the end of inflation, we see that the oscillation condition in Eq. ^ is satisfied just a 
few Hubble times after the end of inflation. Here, we have defined x en d = log(0o/0 e nd)- 



3 The condition that the curvaton keeps following the attractor solution by the end of inflation is given 
by m a < m /(2(log(</> o /<?(> C nd)) 1/2 )- 
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Thus, it is reasonable to assume that the curvaton oscillation starts before the reheating 
of the universe. 

When the curvaton field starts oscillation, the potential of the curvaton is governed 
by the quartic term. In this period, the energy density of the curvaton oscillation scales 
as the energy density of the radiation. The power spectrum of the curvature perturbation 
in such a case is given in Ref. [H] (see also the appendix [A|), 



*<*> = (4+iOT • (23) 

Here, ct OS c and 07- are the field values of the curvaton at the beginning of the curvaton 
oscillation and at the horizon exit of the wave number k during inflation, respectively. As 
we have mentioned above, the curvaton starts oscillation just after the end of inflation, 
and hence, cr OS c ~ Cj, while <7fc is given by <Jk = m^/y/^X log(0o/0fc) ^ 2 with 0^ being 
the field value of the inflaton at the horizon exit of the wave number k. As a result, the 
power spectrum is given by 

— 9 &(t^)(w)(£) 2 (t) 3 ' <*> 

where we have defined Xk = log(0o/0fc) and assumed that r dec "C 1. The nonlinearity 
parameter f^L is also given by 

16(l + r dcc ) 4 + 3r dec ( n 
JNL — 7-: , » r + 



'"dec 

(4 + 3r dec ) r dec \ ct 

4 / -^end 



2 

osc 



(25) 

'"dec V %k 

for r dec <C 1 (see also the appendix[A|). 

Finally, let us estimate the energy fraction of the curvaton at the decay time, r dec . 
The energy densities of the radiation and the curvaton at the curvaton decay time are 
given by 



Pr I dec P<j> I end ^ 



x 3 / 
a cnd \ a R 



a R J \ a dec 



P<r|dec = Palend X ( — ) ( — ^ ) , (26) 

Q>2 J \Ctde 
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where a's are the scale factors of the universe. The meanings of the subscripts are under- 
stood. Therefore, r^ ec can be represented in terms of the scale factors as 



T d e 



^dec 



«2 



^end 



with 



Tend 



Pa 
Pr 



(27) 



(2c 



end 



The first factor is given by ctdec/aR = Tn/Td ec , since the temperature of the universe is 
inversely proportional to the scale factor after the reheating. Between t osc and t 2 ? the 
amplitude of the curvaton field is inversely proportional to the scale factor, and hence, we 
obtain a osc /ci2 = cr 2 /crj = 2m cr /m <} i, x a/ 2 - end- The third factor is, on the other hand, given 
by 



^end 



tt2 
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HI, 
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20 



■^end • 



end 



The curvaton energy fraction at the end of inflation, r en d, is simply given by 

Tend 1 n TTi — , I ~TT I X 



32A V Af, 



'pi 



end 



Putting all of them together, rd ec is given by 

Tdec 



5 m a m^ T R _ 1/2 



12A M 2 T, 



dec 



0.26 x 



"end 
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11 



GeV/ 



T, 



R 



10 12 GeV 



10 6 GeV\ /10 



T, 



dec 



.10 12 GeV 

As a result, by substituting the above rd ec into Eq. (24), we obtain 

. / m \ 2 / Tn 1 \ 2 / Tn \ 2 

P C (k) 



X 



end 



(29) 

(30) 

(31) 

1/2 



X 



10" 



A 



10 n GeV 
7 
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•^end 



10 12 GeV/ V10 12 GeV 
60 



10 6 GeV\ 

Tdec J 



(32) 



while the nonlinearity parameter in Eq. (25) is of O(10) for the above parameter region. 
Therefore, we find that the curvaton model where the initial amplitude of its oscillation is 
set by the attractor behavior can explain the observed power spectrum while predicting a 
rather large nonlinearity parameter. We should again emphasize that the above result is 
insensitive to the initial condition of the curvaton field as long as it follows the attractor 
behavior during inflation. 
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3.3 Cr att (0end) < CT2 

For a larger curvaton mass, m a > m0/(2y / a; end ), the curvaton field reaches a 2 during 
inflation^] In this case, the curvaton stops following the attractor behavior during inflation 
and stays at a ~ ct 2 until it starts oscillation. When the Hubble parameter decreases 
below the curvaton mass (see Eq. ([7])), the curvaton field starts oscillation as a massive 
field. The power spectrum for such a curvaton is given in the appendix[A] (for n = 3, 
V = I Aa 4 + \mla 2 and cr osc = a 2 ), 

r>f,\ r dcc 9A f H kY 3r dec A f 171 A 2 AT 



(4 + 3r dec ) 2 2ml \^ J 4(4 + 3r dcc ) 2 vr 2 Vm. 



?"dec \ 2 / A \ / N k \ ( TTlcj, 



2.5x10- ^ [— ^ ^ . (33) 



2 



0.1 y v i( 3 6 J V 60 y \m c 

The nonlinearity parameter is also given by, 

40(1 + r dec ) 5(4 + 3r dcc ) 86 

JNL ^ ;: , V + 



3r dcc (4 + 3r dcc ) 6r dcc 225 

622 



135 r dec 

By remembering that the energy densities at the curvaton decay time are given by 



(34) 



3 / \ 4 

«end \ (lR 



lend 

a R J V^dec, 

Paldcc - P<j|end(-^-) , (35) 
V a dec / 

we find that the energy ratio can be expressed in terms of the scale factor of the universe, 

(<^osc \ / C-dec \ /o/->\ 
Oend/ V a R J 

As a result, by using the ratios of the scale factors obtained in the previous section, we 
find 

10 /m„\ 2 T R _ /10 6 \ / m„ \ 2 / Tn \ (2 x 10 6 GeV 



ice 



9X\M p J T dec ~ ' X V A J U0 12 GeW Vl0 12 GeV 



R -0ix U^Vl^^M^^) .(37) 



- dec 



4 For m a > 77i0, the curvaton field is settle to its minimum at a = during inflation, and hence, we 
assume m a < in our analysis. 
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Therefore, by putting r^ ec into the above power spectrum in Eq. (33), we obtain 



x m , (38) 



while obtaining a rather large nonlinearity parameter /nl = 0(10) (see Eq. pip). 
3.4 Implications on particle physics models 

Finally, let us discuss some implications of the model. First, let us comment on the decay 
temperature of the curvaton field. As we have discussed in section[2] based on the simplest 
curvaton model, the decay temperature cannot be very high if we require a rather large 
nonlinearity parameter, /tvl = O(10). Here, we show that the similar conclusions are 
reached in the present model. To see such a constraint explicitly, let us express the decay 



temperature Td ec in terms of P^ and f^L by using Eq. (34), (37) and rt38 



p 311N k T R ml | 

- dec 



1296tt 2 m^M^P^f, 



NL 



This shows that the required decay temperature is again rather suppressed for a sizable 
nonlinearity parameter, although it is slightly larger than the one given in section[2] by 
about a factor of 4. 

Secondly, let us discuss a possible interrelation of the present model to another well- 
motivated new physics model, the see-saw mechanism [21] . For that purpose, let us re- 



member that the required quartic coupling A can be expressed by (see Eqs. (32) and (38)) 



97200tt 



96721AV JV ^ ^ v ' \m+ t 

6X W%Y(^| (^)(_*-^| ■ (40) 




Eqs. (32) and (38 1 show that the higher decay temperature is allowed for the higher curvaton mass. 
Since we are interested in the upper bound on the required decay temperature, we concentrate on the 
region m$/(2y/ x cn d) < m a < in the following arguments. 
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If we assume that the mass of the curvaton field is generated by the vacuum expectation 
value of a field X in a supersymmetric model via a superpotential, 

W = gXHZ , a = -^=Re(E) , (41) 

the mass and the quartic coupling of the curvaton are given by m a = g (X) and A = g 2 , 
respectively. Thus, for example, to realize A ~ 10 -6 and m a ~ 10 12 GeV, the vacuum 
expectation value of X is (X) ~ 10 15 GeV. It is quite suggestive that this value is close to 
the mass scale of the seesaw mechanism. Actually, it has been discussed that the right- 
handed scalar neutrinos can play a role of the curvaton field in the supersymmetric seesaw 
mechanism [22] . Detailed studies of the attractor behavior right-handed scalar neutrinos 
as well as the compatibility with leptogenesis will be given elsewhere [23J • 

Thirdly, let us comment on an implication on compensated isocurvature perturbations 
[21]. If the dark matter abundance is created before the curvaton decays, unacceptably 
large matter isocurvature perturbations are generated and contradicts with the recent 
observations of the CMB [6J. This can happen in the models of dark matter such as the 
gravitino produced by thermal scatterings or its decay products [2S], the Q-ball or its decay 
products [261127], the WIMPZILLA [2H] and the primordial black hole [29]. However, the 
constraint can be relaxed if the baryon asymmetry is created as the curvaton decays, 
since the dark matter and the baryon isocurvature perturbations compensate for each 
other. Such scenario is possible if the right-handed scalar neutrino is the curvaton. The 
compensation requires r dcc = Q b /Q DM ~ 0.2 [2"4"] . where fib and f2 DM are the energy density 
fraction of the baryon and the dark matter, respectively, and hence the nonlinearity 
parameter is predicted in each curvaton models. For a curvaton model with a quadratic 
potential, from Eq. ([6]), /nl — 8. In the case of an attractor curvaton with a quartic 
potential, from Eqs. (25) and (34), /nl can be as large as 40. 

As a final remark, let us comment on the Hubble induced mass. In this paper, we have 
assumed that the Hubble induced mass is negligible. In supergravity, scalar fields obtain 
soft masses as large as the Hubble scale if inflation is driven by F-terms [30] , unless there 
exists a tuning of few percents, a Heisenberg symmetry [31] . a discrete R symmetry [32] 
or a shift symmetry [33], or scalar fields are pseudo Nambu-Goldstone bosons [31]. It 
will be interesting to construct curvaton models with an attractor behavior but without 
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fine-tuning in supergravity. 
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A General formulas for the power spectrum and the 
non-Gaussianity 

In this appendix, we summarize the formulas for the curvature perturbation and the 
non-Gaussianity given in Ref. [H]. The assumptions are 

• The energy density of the curvaton oscillation is negligible compared to the total 
energy density at least until the onset of the curvaton oscillation. 

• The curvaton field starts oscillation suddenly at t = t osc and its energy density scales 
by a" n in terms of the scale factor of the universe a. 

■ At some point before the curvaton decays, the quadratic term dominates and the 
energy density of the curvaton begins to scale by a -3 . 

• The curvaton oscillation decays suddenly and the curvature perturbation is fixed at 
that point. 

• The curvature perturbation from the fluctuation of the inflaton field is negligible. 

• The energy density of the inflaton field scales by a~ 3 after the end of inflation until 
the reheating. 

• The reheating completes instantaneously at t = £r before the curvaton decays. 
With these assumptions, the power spectrum of the curvature perturbation is given by 



7> c (fc)* 



4 + 3r dec 



'"dec 



(1 - X{a osc )) 



V^Vosc) p V^Vosc) _ 3X(a osc ) \ Ih 
V'(a k ) \n V(a osc ) a osc J 2ir 



(42) 
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where the function X is defined by 



X{a) 



1 



aV"(a) 



- 1 



2(c-3) V V'(a) 

5 (t R < t osc ) , 
9/2 (t osc < . 



(43) 



The nonlinearity parameter is given by 



f NL {k) 



40(l + r dec ) , 5(4 + 3r dec ) {3V'(a. 



3r dec (4 + 3r dec ) 



+ 



X'(<T n 



1 - X(cr osc ) 
3 {V'(a nRC 



3 V'(<r n 



n V(a, 
3X(a, 



n \ V(a, 
V"(a c 




-i 



x 



n V{a, 



V'(a c 



- [I- X{a osc )) 

V\(T osc ) 



(44) 
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